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Abstract 

The ATLAS experiment at the LHC has measured the production cross section of events with two 
isolated photons in the final state, in proton-proton collisions at y/s = 7 TeV. The full data set collected 
in 201 1, corresponding to an integrated luminosity of 4.9 fb~\ is used. The amount of background, 
from hadronic jets and isolated electrons, is estimated with data-driven techniques and subtracted. 
The total cross section, for two isolated photons with transverse energies above 25 GeV and 22 GeV 
respectively, in the acceptance of the electromagnetic calorimeter (|r?| < 1.37 and 1.52 < \r]\ < 2.37) 
and with an angular separation AR > 0.4, is 44.0^4 2 pb. The differential cross sections as a function 
of the di-photon invariant mass, transverse momentum, azimuthal separation, and cosine of the polar 
angle of the largest transverse energy photon in the Collins-Soper di-photon rest frame are also 
measured. The results are compared to the prediction of leading-order parton-shower and next-to- 
leading-order and next-to-next-to-leading-order parton-level generators. 
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ABSTRACT: The ATLAS experiment at the LHC has measured the production cross section 
of events with two isolated photons in the final state, in proton-proton collisions at y/s = 
7 TeV. The full data set collected in 2011, corresponding to an integrated luminosity of 
4.9 fb" 1 , is used. The amount of background, from hadronic jets and isolated electrons, 
is estimated with data-driven techniques and subtracted. The total cross section, for two 
isolated photons with transverse energies above 25 GeV and 22 GeV respectively, in the 
acceptance of the electromagnetic calorimeter (|r/| < 1.37 and 1.52 < \rj\ < 2.37) and with 
an angular separation AR > 0.4, is 44. pb. The differential cross sections as a function 
of the di-photon invariant mass, transverse momentum, azimuthal separation, and cosine of 
the polar angle of the largest transverse energy photon in the Collins-Soper di-photon rest 
frame are also measured. The results are compared to the prediction of leading-order parton- 
shower and next-to-leading-order and next-to-next-to-leading-order parton-level generators. 
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1 Introduction 

The measurement at the LHC of the production cross section, in pp collisions, of two 
isolated photons not originating from hadronic decays, pp — > 77 + X, provides a tool 
to probe perturbative Quantum Chromo dynamics (QCD) predictions and to understand 
the irreducible background to new physics processes involving photons in the final state. 
These processes include Higgs boson decays to photon pairs (H — > 77) or graviton decays 
predicted in some Universal Extra-Dimension models [1, 2]. 

Recent cross section measurements for di-photon production at hadron colliders were 
performed by the D0 [3] and CDF [4] collaborations at the y/s = 1.96 TeV Tevatron pp 
collider, and by ATLAS [5] and CMS [6] using ^/s = 7 TeV pp collisions recorded at the 
LHC in 2010. 

In this paper, the production cross section of two isolated photons with transverse 
energies (E^) above 25 GeV and 22 GeV respectively, in the acceptance of the ATLAS elec- 
tromagnetic calorimeter (\rj\ < 1.37 and 1.52 < | -^7 1 < 2.37) and with an angular separation 
AR > 0.4, is measured. The results are obtained using the data collected by the ATLAS 
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experiment in 2011, which corresponds to an integrated luminosity 1 of (4.9 ± 0.2) fb , 
thus increasing the sample size by more than a factor of 100 compared to the previous 
measurement. The transverse energy thresholds for the two photons are higher than in the 
previous measurement (16 GeV). 

The integrated di-photon production cross section is measured, as well as the differen- 
tial cross sections as a function of four kinematic variables: the di-photon invariant mass 
(m 77 ), the di-photon transverse momentum (pT,yy), the azimuthal 2 separation between the 
photons in the laboratory frame (A^~~), and the cosine of the polar angle of the highest 
photon in the Collins-Soper di-photon rest frame (cos#* 7 ) [9]. The first distribution is of 
obvious interest for resonance searches; the second and the third provide important infor- 
mation in the study of higher-order QCD perturbative effects and fragmentation, especially 
in some specific regions such as the small A0 77 limit; the fourth can be used to investi- 
gate the spin of di-photon resonances. The results are compared to the predictions from: 
parton-shower Monte Carlo generators, PYTHIA [10] and SHERPA [11]; parton-level calcu- 
lations with next-to-leading-order (NLO) QCD corrections using the DlPHOX [12] program 
complemented by GAMMA2MC [13]; and at next-to-next-to-leading-order (NNLO), using 
27NNLO [14]. The contribution from the di-photon decays of the particle recently discov- 
ered by ATLAS [15] and CMS [16] in the search for the Standard Model Higgs boson is 
not included in the theoretical calculations. It is expected to contribute around 1% of the 
signal in the 120< m 77 <130 GeV interval, and negligibly elsewhere. 

2 The ATLAS detector 

ATLAS [17] is a multipurpose detector with a forward-backward symmetric cylindrical 
geometry and nearly 4tt coverage in solid angle. The most relevant subdetectors for the 
present analysis are the inner tracking detector (ID) and the calorimeters. 

The ID consists of a silicon pixel detector and a silicon microstrip detector covering 
the pseudorapidity range \r]\ < 2.5, and a straw tube transition radiation tracker covering 
\r/\ < 2.0. It is immersed in a 2 T magnetic field provided by a superconducting solenoid. 
The ID allows efficient reconstruction of converted photons if the conversion occurs at a 
radius of up to ~ 0.80 m. 

The electromagnetic calorimeter (ECAL) is a lead/liquid-argon (LAr) sampling calo- 
rimeter providing coverage for \rj\ < 3.2. It consists of a barrel section (\r]\ < 1.475) 
and two end-caps (1.375 < \r]\ < 3.2). The central region (\rj\ < 2.5) is segmented into 
three longitudinal layers. The first (inner) layer, covering \rj\ < 1.4 in the barrel and 
1.5 < \r}\ < 2.4 in the end-caps, has high granularity in the rj direction (between 0.003 and 

The 3.9% uncertainty in the integrated luminosity for the complete 2011 data set is based on the 
calibration described in refs. [7, 8] including an additional uncertainty for the extrapolation to the later 
data-taking period with higher instantaneous luminosity. 

2 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in 
the centre of the detector and the 2-axis coinciding with the axis of the beam pipe. The i-axis points from 
the IP to the centre of the LHC ring, and the y-axis points upward. Cylindrical coordinates (r,(j>) axe used 
in the transverse plane, <j> being the azimuthal angle around the beam pipe. The pseudorapidity is defined 
in terms of the polar angle 8 as rj = — lntan(6 l /2). 
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0.006 depending on r/), sufficient to provide event-by-event discrimination between single- 
photon showers and two overlapping showers from a 7r° decay. The second layer, which 
collects most of the energy deposited in the calorimeter by the photon shower, has a cell 
granularity of 0.025x0.025 in rj x <j). The third layer is used to correct high energy showers 
for leakage beyond the ECAL. In front of the electromagnetic calorimeter a thin presampler 
layer, covering the pseudorapidity interval \r)\ < 1.8, is used to correct for energy loss before 
the ECAL. 

The hadronic calorimeter (HCAL), surrounding the ECAL, consists of an iron/scin- 
tillator tile calorimeter in the range \r]\ < 1.7, and two copper/LAr calorimeters spanning 
1.5 < |?7| < 3.2. The ECAL and HCAL acceptance is extended by two LAr forward 
calorimeters (using copper and tungsten as absorbers) up to \r]\ < 4.9. 

Di-photon events are recorded using a three-level trigger system. The first level, im- 
plemented in hardware, is based on towers defined with a coarser granularity (0.1 x 0.1 in 
7] x (/>) than that of the ECAL. They are used to search for electromagnetic deposits in 
7] x <j) regions of 2 x 1 and 1x2 towers, within a fixed window of size 2x2 and with a 
transverse energy above a programmable threshold. The second- and third-level triggers 
are implemented in software and exploit the full granularity and energy calibration of the 
calorimeter to refine the first-level trigger selection. 



3 Data and Monte Carlo samples 

The data set analysed consists of the 7 TeV proton-proton collisions recorded by the ATLAS 
detector in 2011. Only events where the beam conditions are stable and the trigger system, 
the tracking devices, and the calorimeters are operational, are considered. 

Monte Carlo (MC) samples are produced using various generators as described below. 
Particle interactions with the detector material and the detector response are simulated 
with GEANT4 [18]. The events are reconstructed with the same algorithms used for collision 
data. More details of the event generation and simulation infrastructure are provided in 
ref. [19]. 

Simulated di-photon events are generated with both PYTHIA 6.4.21 and SHERPA 1.3.1. 
PYTHIA uses the modified leading-order MRST2007 [20] parton distribution functions 
(PDFs) while SHERPA uses the CTEQ6L1 [21] PDFs. The event-generator parameters 
are set according to the ATLAS AMBT2 [22] tune. Photons originating from the hard 
scattering and quark bremsstrahlung are included in the analysis. The MC di-photon sig- 
nal is generated with a photon threshold of 20 GeV; one million events are produced 
both with PYTHIA and SHERPA. They are used to model the transverse isolation energy 
(see section 4) distribution of signal photons, to compute the reconstruction efficiency and 
to study the systematic uncertainties on the reconstructed quantities. Background 7-jet 
events are generated using Alpgen [23] with the CTEQ6L1 PDF set. 
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4 Event selection 



Events are collected using a di-photon trigger with a nominal transverse energy threshold 
of 20 GeV for both photon candidates. The photon trigger objects are required to pass a 
selection based on shower shape variables computed from the energy deposits in the second 
layer of the electromagnetic calorimeter and in the hadronic calorimeter. The requirements 
are looser than the photon identification criteria applied in the offline selection. In order 
to reduce non-collision backgrounds, events are required to have a reconstructed primary 
vertex with at least three associated tracks and consistent with the average beam spot 
position. The signal inefficiency of this requirement is negligible. 

Photons are reconstructed from electromagnetic energy clusters in the calorimeter and 
tracking information provided by the ID as described in ref . [24] . Photons reconstructed near 
regions of the calorimeter affected by read-out or high-voltage failures are not considered. 3 
The cluster energies are corrected using an in-situ calibration based on the Z boson mass 
peak [25], and the determination of the pseudorapidities is optimized using the technique 
described in ref. [15]. In order to benefit from the fine segmentation of the first layer of 
the electromagnetic calorimeter to discriminate between genuine prompt photons and fake 
photons within jets, the photon candidate pseudorapidity must satisfy \r\\ < 1.37 or 1.52 < 
|?7| < 2.37. We retain photon candidates passing loose identification requirements, based on 
the same shower shape variables - computed with better granularity and resolution - and 
the same thresholds used at trigger level. The highest-Ii/p ("leading") and second highest- 
-Et ("subleading") photons within the acceptance and satisfying the loose identification 
criteria are required to have £7t,i > 25 GeV and £r,2 > 22 GeV, respectively. The angular 
separation between the two photons, AR = sj (Ar/) 2 + (A(f>) 2 , is required to be larger than 
0.4, in order to avoid one photon candidate depositing significant energy in the isolation 
cone of the other, as defined below. 

Two further criteria are used to define the signal and background control regions. Firstly 
the tight photon selection [24] (abbreviated as T in the following) is designed to reject 
hadronic jet background, by imposing requirements on nine discriminating variables com- 
puted from the energy leaking into the HCAL and the lateral and longitudinal shower 
development in the ECAL. Secondly the transverse isolation energy E^° is computed from 
the sum of the positive-energy topological clusters with reconstructed barycentres inside 
a cone of radius AR = 0.4 around the photon candidate. The algorithm for constructing 
topological clusters suppresses noise by keeping only those cells with a significant energy 
deposit and their neighbouring cells. The cells within 0.125 x 0.175 in rj X (f) around the 
photon are excluded from the calculation of E^°. The mean value of the small leakage of 
the photon energy from this region into the isolation cone, evaluated as a function of the 
photon transverse energy, is subtracted from the measured value of E^° (meaning that 
can be negative). The typical size of this correction is a few percent of the pho- 
ton transverse energy. The measured value of E^° is further corrected by subtracting the 
estimated contributions from the underlying event and additional pp interactions. This cor- 

3 This requirement leads to a typical loss of 0.8% to 1.4% on the photon reconstruction efficiency, de- 
pending on the data-taking period. 
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rection is computed on an event-by-event basis, by calculating the transverse energy density 
from low-transverse-momentum jets, as suggested in refs. [26, 27]. The median transverse 
energy densities of the jets in two r\ regions, \r)\ < 1.5 and 1.5 < \rj\ < 3.0, are computed 
separately, and the one for the region containing the photon candidate pseudorapidity is 
multiplied by the total area of all topological clusters used in the calculation of the isola- 
tion variable in order to estimate the correction. Signal photons are required to pass the 
tight selection ("tight photons") and the isolation requirement I, —4 < < 4 GeV. A 
total of 165 767 pairs of tight, isolated photons are selected. The non-tight (T) photon 
candidates are defined as those failing the tight criteria for at least one of the shower-shape 
variables that are computed from the energy deposits in a few cells of the first layer of the 
electromagnetic calorimeter adjacent to the cluster barycentre. Photon candidates with 

4 < < 8 GeV are considered non-isolated (I). 

5 Signal yield extraction 

After the selection, the main background is due primarily to 7-jet and secondarily to di-jet 
(jj) final states, collectively called "jet background" in the following. Two methods, the 
two-dimensional sidebands and the two-dimensional fit, already exploited in ref. [5], are 
used to perform an in-situ statistical subtraction of the jet background from the selected 
photon candidate pairs, as described in section 5.1. 

After the jet background contribution is subtracted, a small residual background con- 
tamination arises from events where isolated electrons are misidentified as photons. This 
contribution is estimated as described in section 5.2. 

5.1 Jet background subtraction 

Both the two-dimensional sidebands and the two-dimensional fit methods use the photon 
transverse isolation energy and the tight identification criteria to discriminate prompt pho- 
tons from jets. They rely on the fact that the correlations between the isolation and the 
tight criteria in background events are small, and that the signal contamination in the 
non-tight or non-isolated control regions is low. 

The two-dimensional sidebands method counts the numbers of photon candidate pairs 
where each of the candidates passes or fails the tight and the isolation criteria. Four 
categories are defined for each photon, resulting in 16 categories of events. The inputs 
to the method are the numbers of events in the categories and the signal efficiencies of 
the tight and isolation requirements. The correlation between these two requirements is 
assumed to be negligible for background events. The method allows the simultaneous 
extraction of the numbers of true di-photon signal, 7J, j7 4 and jj background events, and 
the tight and isolation efficiencies for fake photon candidates from jets ("fake rates"). The 
expected number of events in each category is written as a function of the parameters 
(yields, efficiencies, fake rates and correlation factors) and the system of 16 equations is 
solved with a \ 2 minimization procedure. This method is an extension of the one used in 

4 Here and in the following, 7j Q7) denotes the events where the leading (subleading) candidate is a true 
photon, and the other candidate a true jet. 
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our previous di-photon analysis [5]. It allows the extraction of different isolation fake rates 
for jets in J7 or jj events as well as a correlation factor for the isolation of jet pairs. 

The two-dimensional fit method consists of an extended maximum likelihood template 
fit to the two-dimensional distribution of the transverse isolation energies E!^\ and E^° 2 
of the two photon candidates in events belonging to the T-T sample, i.e. where both 
photons satisfy the tight identification criteria. The fit is performed in the isolation range 
—4 < < 8 GeV (i = 1, 2). The correlations between the transverse isolation energies of 
the two candidates in di-photon, 7J, and J7 events are found to be negligible in MC samples, 
and the products of two one-dimensional templates for E^ and E^° 2 are used for each of 
the three event species. For the jj component, large correlations are observed in data, and 
a two-dimensional template is used. The two-dimensional fit is described in detail in our 
previous paper [5]. There are two differences between the present and previous analyses: 
the use now of binned distributions instead of smooth parametric functions for the photon 
and jet templates, and the correction for signal leakage in the background templates, as 
described below. 

The transverse isolation energy distributions of the signal photons and the correspond- 
ing efficiencies of the signal requirement —4 < E™° < 4 GeV are obtained from the SHERPA 
di-photon sample, separately for the leading and the subleading candidates. In the two- 
dimensional fit method, the templates are shifted by +160 and +120 MeV respectively in 
order to maximize the likelihood, as determined from a scan as a function of the shifts. 
These values are also used to compute the signal efficiencies of the isolation requirement 
needed in the two-dimensional sidebands method. Shifts of similar size between ATLAS 
data and MC simulation have been observed in the transverse isolation energy distribution, 
computed with the same technique (based on topological clusters inside a cone of radius 
0.4), of electron control samples selected from Z — > ee decays with a tag-and-probe method. 
The E^° distributions of prompt photons in 7j and j7 events are assumed to be identical to 
that of prompt photons in di-photon events, as found in simulated samples. The tight iden- 
tification efficiencies for prompt photons, needed in the two-dimensional sidebands method 
and in the final cross section measurement, are estimated using the same di-photon MC 
sample. The shower shape variables are corrected for the observed differences between data 
and simulation in photon-enriched control samples. Residual differences between the ef- 
ficiencies in the simulation and in data are corrected using scale factors determined from 
control samples of photons from radiative Z boson decays, electrons selected with a tag- 
and-probe technique from Z — > ee decays, and photon-enriched control samples of known 
photon purity [28]. After applying these corrections, the photon identification efficiency in 
the simulation is estimated to reproduce the efficiency in data to within 2%. For the two- 
dimensional fit, the transverse isolation energy template of the leading (subleading) jet in j7 
(7j) events is extracted directly from data where one candidate passes the non-tight and the 
other passes both the tight identification and isolation (TI) requirements. For jj events, the 
two-dimensional template is obtained from data in which the two candidates are required 
to be non-tight. The correlation is found to be about 8%. The jet background templates 
are corrected for signal leakage in the control samples, estimated from the SHERPA sample. 

Figure 1 shows the projections of the two-dimensional fit to the transverse isolation 
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energies of the leading and subleading photon candidates. The yields for each of the four 
components extracted with the two-dimensional sidebands method and the two-dimensional 
fit are given in table 1. The di-photon purity is around 68% and the di-photon yields agree 
within 1.5% between the two methods. 




Figure 1. Projections of the two-dimensional fit to the transverse isolation energies of the two 
photon candidates: leading photon (left) and sub- leading photon (right). The photon templates 
from Sherpa are shifted by +160 MeV (+120 MeV) for the leading (subleading) photon. Solid circles 
represent the observed data. The (black) solid line is the fit result, the (violet) dash-dotted curve 
shows the 77 component. The (red) dotted line shows in the left (right) figure the contribution from 
7J O7) events. In both figures, the (blue) dashed line represents a broad background component 
in the photon candidates' sample: for the leading candidate this is due to j7 and jj final states, 
whereas for the sub-leading candidate it comes from 7] and jj final states. 
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Table 1. Total yields for two candidates satisfying the tight identification and the isolation 
requirement —4 < < 4 GeV. Both statistical and total systematic uncertainties are listed. 

To obtain the differential signal yields as a function of the di-photon kinematic vari- 
ables, such as m TT , Pt,77> A0 77 and cos#* 7 , the above methods are applied in each bin of 
the variables. Figure 2 shows the differential spectra of the signal and background com- 
ponents obtained with the two-dimensional fit. In some regions of the di-photon spectra, 
discrepancies with the two-dimensional sidebands results are larger than those observed 
for the integrated yield. The results from the two-dimensional fit are used to extract the 
nominal cross sections, while differences between the results obtained with the two methods 
are included in the final systematic uncertainty. 

Several sources of systematic uncertainty on the signal yield, estimated after the jet 
background subtraction, are considered. The dominant uncertainty originates from the 
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Figure 2. Differential spectra in data (solid circles) and from the two-dimensional fit, for the 
77 (hollow histogram), 7j+j7 (light solid histogram), and jj (dark solid histogram) contributions. 
The spectra are shown for the following di-photon variables: m 77 (top left), PT.77 (top right), 
A(/> 77 (bottom left), cos 6* (bottom right). 



choice of the background control regions and accounts for both the uncertainty on the 
background transverse isolation energy distribution and its correlation with the identifica- 
tion criteria. It is first estimated by varying the number of relaxed criteria in the non-tight 
definition. For the integrated di-photon yield, the effect is found to be ^g%. In some bins 
of the ?7i 77 and Pt,^ spectra where the size of the control samples is small, neighbouring 
bins are grouped together to extract the jet background templates. Since the background 
transverse isolation energies depend mildly on these kinematic variables, a systematic un- 
certainty is evaluated by repeating the yield extraction with jet templates from the adjacent 
groups of neighbouring bins. The uncertainty on the estimated signal yield is at most ±9%. 

In the nominal result, the photon isolation templates are taken from the SHERPA di- 
photon sample. A systematic uncertainty is evaluated by using alternative templates from 
the PYTHIA di-photon sample, and from data. The data-driven template for the leading 
(subleading) photon is obtained by selecting events where the requirement Ei^° < 8 GeV is 
removed for the leading (subleading) photon candidate, and normalizing the leading (sub- 
leading) photon isolation distribution in T-TI (TI-T) events, where the leading (sublead- 
ing) candidate fails the tight identification while the other candidate passes tight identifica- 
tion and isolation criteria, to the isolation distribution of leading (subleading) candidates in 
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T-T events in the 7 < E^° < 17 GeV region. The difference between the two distributions 
is used as an estimate of the photon distribution. The PYTHIA di-photon sample exhibits 
higher tails than SHERPA at large values of E™°. The data-driven template, on the other 
hand, is characterized by smaller tails than the SHERPA template, since it is obtained by 
assuming that the isolation region above 7 GeV is fully populated by background. The 
corresponding uncertainty on the signal yield is estimated to be of the integrated di- 
photon yield. It is rather uniform as a function of m 77 , PT,-yy, A^ 77 and cos#* 7 and always 
below 4%, except at very low m 77 where it reaches ±5%. The photon isolation template 
is, to a large extent, independent of the variables under study. Repeating the background 
subtraction procedure using photon isolation templates extracted in bins of the di-photon 
variable under study leads to variations of the estimated signal yield within ^4%. 

Other systematic effects have been considered, and found to be smaller than those 
previously discussed. The bias created by neglecting the dependence of the identification 
and isolation efficiencies on r\ and E^ is estimated to be of +0.02% and -0.3% respectively. 
The effect of assuming identical templates for photons in di-photon and in 7-jet events is 
evaluated by using instead templates from ALPGEN 7-jet samples for photons in the 7j 
and j7 components. The uncertainty on the shifts applied to the MC photon templates 
(±10 MeV for the leading photons and ±5 MeV for the subleading ones, as determined 
from the scan) is propagated to the di-photon yields. The impact of the identification 
efficiencies on the signal leakage correction is estimated by neglecting in the simulation 
the correction factors nominally applied to the shower shape variables to account for the 
observed differences between data and MC simulation. These effects produce systematic 
uncertainties of at most 0.5% on the differential spectra. Finally, no significant effect is 
observed due to the imperfect modelling of the material in front of the calorimeters. 



5.2 Electron background subtraction 

Isolated electrons from W or Z boson decays can be misidentified as photons, since the 
two particles (e and 7) generate similar electromagnetic showers in the ECAL. Usually a 
track is reconstructed in the inner detector pointing to the electron ECAL cluster, thus 
isolated electrons misidentified as photons are mostly classified as converted candidates. 
Pairs of misidentified, isolated electrons and positrons (ee) from processes such as Drell- 
Yan, Z —7- ee, WW — > eveis, or of photons and e ± from diboson production (-yW — > jeu, 
"fZ —7- 7ee), provide a background that cannot be distinguished from the di-photon signal 
based on the photon identification and isolation variables and must therefore be estimated 
in a separate way. The same procedure exploited in ref. [5], based on the number of 7e 
(A^ye) and ee (N ee ) events observed in data, is used to estimate their contributions to the 
di-photon yield iV 77 after jet background subtraction. 

For a given bin i of the variable X {X = m 77 , Pt,77, A^ 77 , or cos#* 7 ), the signal 
component Nyf in the A 77 sample can be evaluated: 

jusig _ ^77 ~ [fe-^y^ye ~ (/e-^) 2 Nee] , , 

iV 77 ~ 7] f f ^2 

K 1 Je-^^J^-^e) 
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The fake rates / e ->-y and / 7 _j. e are measured using Z boson decays in data. Z — > ee de- 
cays are used to estimate / e ^ 7 as N^ e /(2N^ e ), where N^ e and N^ e are the numbers of 
7e and ee pairs with invariant mass within 1.5 a of the Z boson mass. Z — > "fee decays 
are similarly used to estimate / 7 -j. e as N^ ee /N^ ee . The numbers of continuum background 
events are estimated from the sidebands of the ee, je, eee or 7ee invariant mass distribu- 
tions (51 - 61GeV and 121 - 131 GeV), and subtracted from N^ e , N^ e , N^ ee and N^ ee , 
respectively. Electrons must satisfy identification criteria based on their shower shape in 
the electromagnetic calorimeter, quality criteria for the associated track in the ID, and 
an isolation requirement < 4 GeV. The measured fake rates, including statistical and 
systematic uncertainties, are / e _> 7 = 0.062^q'q^q and / 7 _>. e = 0.038^q'qq7, where the system- 
atic uncertainty is dominated by the dependence on the transverse energy of the candidate 
photon. Other sources include the uncertainties on N^ e , N? e , N^ ee and N^ ee which are 
evaluated by changing the definition of the Z boson mass window to ±2cx and ±1<7, and 
shifting the sidebands by ±5 GeV. The fraction of electron background as a function of 
m 77 , j*t,77) A</Vp an d cos#* 7 is shown in figure 3. The enhancements at m 77 ~ mz, low 
Pt,77 an d A(/> 77 ~ 7r are due to the large Z boson production cross section. 
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6 Cross section measurement 

This section describes the extraction of the final cross sections. The background-subtracted 
differential spectra are first unfolded to the generated-particle level, to take into account 
reconstruction and selection efficiencies estimated from the simulation, and then divided by 
the integrated luminosity of the data sample and the trigger efficiency relative to the offline 
selection. 

6.1 Efficiency and unfolding 

The background-subtracted differential distributions obtained from the data are unfolded to 
obtain the particle-level spectra by dividing the signal yield in each bin of the di-photon ob- 
servable under study by a "bin-by-bin" correction, which accounts for signal reconstruction 
and selection efficiencies and for finite resolution effects. The bin-by-bin nominal corrections 
are evaluated from the SHERPA di-photon simulated sample as the number of simulated di- 
photon events satisfying the selection criteria (excluding the trigger requirement) and for 
which the reconstructed value of the variable X under consideration is in bin i, divided 
by the number of simulated di-photon events satisfying the nominal acceptance criteria at 
generator- level and for which the generated value of X is in the same bin i. The generator- 
level photon transverse isolation energy is computed from the true four-momenta of the 
generated particles (excluding muons and neutrinos) inside a cone of radius 0.4 around the 
photon direction. The pileup contribution is removed using an analogous method to the 
one for the experimental isolation variable, by subtracting the product of the area of the 
isolation cone and the median transverse energy density of the low-transverse-momentum 
truth-particle jets. 

Alternative corrections are calculated with the PYTHIA di-photon sample or using a 
simulated di-photon sample which contains additional material upstream of the calorime- 
ter. The variations induced on the measured cross sections by the alternative corrections 
are taken as systematic uncertainties, due to the uncertainty on the generated kinematic 
distributions, on the relative fraction of direct and fragmentation di-photon production, 
and on the amount of material in the ATLAS detector. The effect on the total cross section 
is within +5% for m 77 , ±3% for pt,^, -4% for A^ 77 and +3% for cos#* 7 . 

The effect of the uncertainty on the efficiency of the photon identification criteria is 
estimated by varying the identification efficiency in the simulation by its uncertainty [28]. 
The uncertainties on the electromagnetic (photon) energy scale and resolution are also 
propagated to the final measurement by varying them within their uncertainties [25]. The 
effect on the differential cross section is typically ~^_\%. Other uncertainties, related to the 
dependence on the average number of pile-up interactions of the efficiencies of the photon 
identification and transverse isolation energy requirements and to the observed data-MC 
shift in the photon transverse isolation energy distributions, are found to be negligible. 

A closure test has been performed by unfolding the differential spectra of di-photon 
events selected in the PYTHIA signal sample with the bin-by-bin coefficients determined 
using the SHERPA sample, and comparing the unfolded spectra to the truth-level spectra in 
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the same PYTHIA sample. Non-closure effects of at most 2% have been found and included 
in the final systematic uncertainty. 

More sophisticated unfolding methods which account for migrations between bins, ei- 
ther based on the repeated (iterative) application of Bayes' theorem [29] or on a least-square 
minimization followed by a regularization of the resulting spectra [30] have also been in- 
vestigated. The differences between the unfolded spectra obtained with these methods and 
the spectra extracted with the bin-by-bin corrections are negligible compared to the other 
uncertainties and therefore the bin- by-bin method was chosen for the final results. 

6.2 Trigger efficiency correction 

The unfolded spectra are then corrected for the event-level trigger efficiency, defined as 
the fraction of di-photon events - satisfying all the selection criteria - that pass the di- 
photon trigger used to collect the data. The trigger efficiency is measured in data using 
a bootstrap technique from samples selected with fully efficient unbiased triggers with a 
lower threshold, and taking into account kinematic correlations between the photon candi- 
dates. The differences between the measured di-photon trigger efficiency and the efficiency 
estimated with simulated di-photon samples, or with an alternative data-driven technique 
that neglects the correlations between the two photon candidates, have been assigned as 
systematic uncertainties. The total trigger efficiency is then: 

e trig = (97.8l?|(stat.) ± 0.8(syst.)) % (6.1) 

Using di-photon simulated samples, the trigger efficiency has been estimated to be con- 
stant, within the total uncertainties, as a function of the four di-photon observables under 
investigation. 

6.3 Results 

The differential cross sections as a function of m 77 , pr,77> Ac/> 77 , an d cos#* 7 are extracted 
following the unfolding procedure described in section 6.1 and using the trigger efficiency 
quoted in eq. (6.1). The numerical results are listed in Appendix A. 

The integrated cross section is measured by dividing the global 77 yield (obtained after 
subtracting the electron contribution from the two-dimensional fit result in table 1) by the 
product of the average event selection efficiency (from the simulation), trigger efficiency and 
integrated luminosity. The selection efficiency is defined as the number of reconstructed 
simulated di-photon events satisfying the detector-level selection criteria divided by the 
number of generated events satisfying the equivalent truth-level criteria, thus correcting for 
reconstructed events with true photons failing the acceptance cuts. It is computed from 
simulated di-photon events, reweighting the spectrum of one of the four di-photon variables 
under study in order to match the differential background-subtracted di-photon spectrum 
observed in data. Choosing different variables for the reweighting of the simulated events 
leads to slightly different but consistent efficiencies, with an average value of 49.6% and 
an RMS of 0.2%. Including systematic uncertainties on the photon reconstruction and 
identification efficiencies, from the same sources described in section 6.1, the event selec- 
tion efficiency is estimated to be 49. 6^7%. The dominant contributions to the efficiency 
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uncertainty are from the photon identification efficiency uncertainty (±1.2%), the energy 
scale uncertainty (lj'5%), and the choice of the MC generator and the detector simula- 
tion (±0.9%). Negligible uncertainties are found to arise from the energy resolution, the 
isolation requirement (evaluated by shifting the isolation variable by the observed data- 
MC difference) and from the different pile-up dependence of the efficiency in data and MC 
simulation. With an integrated luminosity of (4.9 ± 0.2) fb _1 at y/s = 7 TeV, we obtain 
an integrated cross section of 44.0^42 pb, where the dominant uncertainties are the event 
selection efficiency and the jet subtraction systematic uncertainties. As a cross-check, the 
integrals of the one-dimensional differential cross sections are also computed. They are 
consistent with the measured integrated cross section quoted above. 

7 Comparison with theoretical predictions 

The results are compared both to fixed-order NLO and NNLO calculations, obtained with 
parton-level MC generators (DIPHOX+GAMMA2MC and 27NNLO), and to the generated- 
particle-level di-photon spectra predicted by leading-order (LO) parton-shower MC gener- 
ators used in the ATLAS full simulation (PYTHIA and Sherpa). The contribution from 
the particle recently discovered by ATLAS and CMS in the search for the Standard Model 
Higgs boson is not included in the theoretical calculations: it is expected to be around 
1% of the signal in the 120< m 77 <130 GeV interval, and negligible elsewhere. The main 
differences between the four predictions are the following: 

• 27NNLO provides a NNLO calculation of the direct part of the di-photon production 
cross section, but neglects completely the contribution from the fragmentation com- 
ponent, where one or both photons are produced in the soft collinear fragmentation 
of coloured partons. 

• DlPHOX provides a NLO calculation of both the direct and the fragmentation parts of 
the di-photon production cross section. It also includes the contribution from the box 
diagram ( gg — > 77) , which is in principle a term of the NNLO expansion in the strong 
coupling constant a s , but - due to the large gluon luminosity at the LHC [31] - gives 
a contribution comparable to that of the LO terms. For these reasons, higher-order 
contributions to the box diagrams, technically at NNNLO but of size similar to that 
of NLO terms, are also included in our calculation by using GAMMA2MC. 

• PYTHIA provides LO matrix elements for di-photon production and models the higher- 
order terms through 7-jet and di-jet production in combination with initial-state 
and/or final-state radiation. It also features parton showering and an underlying 
event model; 

• SHERPA has features similar to those of PYTHIA, and in addition includes the di- 
photon higher-order real-emission matrix elements. For this study, up to two addi- 
tional QCD partons are generated. 

The nominal factorization (/xp), renormalization (jUr), and - in the case of DlPHOX 
and GAMMA2MC - fragmentation (fif) scales are set in all cases to the di-photon invariant 
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mass, m 77 . Different PDF sets are used by each program: CT10 NLO [32] for DlPHOX and 
GAMMA2MC, MSTW2008 NNLO [33] for 2 7 NNLO, CTEQ6L1 for Sherpa and MRST2007 
LO* for PYTHIA. The theoretical uncertainty error bands for PYTHIA and SHERPA include 
only statistical uncertainties. The theory uncertainty error bands for the NLO and NNLO 
predictions include in addition PDF and scale uncertainties. PDF uncertainties are esti- 
mated by varying each of the eigenvalues of the PDFs by ±la and summing in quadra- 
ture separately positive and negative variations of the cross section. For DlPHOX and 
GAMMA2MC, scale uncertainties are evaluated by varying each scale to m 77 /2 and 2m 77 , 
and the envelope of all variations is taken as a systematic error; the final uncertainty is 
dominated by the configurations in which the scales are varied incoherently. For 27NNLO, 
the scale uncertainty is evaluated by considering the variation of the predicted cross sections 
in the two cases = m 77 /2, [ip = 2m 77 and [ir = 2m 77 , \ip = m~ n /2. 

Fixed-order predictions calculated at parton level do not include underlying event, pile- 
up or hadronization effects. While the ambient-energy density corrections to the photon 
isolation are expected to remove most of these effects from the photon isolation energy, it is 
not guaranteed that they correct the experimental isolation back to exactly the parton-level 
isolation computed from the elementary-process partons. To estimate these residual effects, 
PYTHIA and SHERPA di-photon samples are used to evaluate the ratio of generator-level 
cross sections with and without hadronization and the underlying event, and subsequently, 
the parton-level cross sections are multiplied bin-by-bin by this ratio. The central value of 
the envelope of the PYTHIA and SHERPA distributions is taken as the nominal correction 
and half of the difference between PYTHIA and SHERPA as the systematic uncertainty. The 
typical correction factor is around 0.95. 

Both PYTHIA and SHERPA are expected to underestimate the total cross section, be- 
cause of the missing NLO (and higher-order) contributions. At low and for A</> 77 near 
7T where multiple soft gluon emission is important, PYTHIA and SHERPA are expected to 
better describe the shape of the differential distributions, thanks to the effective all-order 
resummation of the leading logs performed by the parton shower. On the other hand, in the 
same regions fixed-order calculations are expected to exhibit infrared divergences. Finally, 
27NNLO is expected to underestimate the data in regions populated by the contribution 
from fragmentation (low A<^> 77 and m 77 , and cos#* 7 ~l). 

The total cross section estimated by PYTHIA and SHERPA with the ATLAS simu- 
lation settings is 36 pb, and underestimates the measured cross section by 20%. The 
DIPHOX+GAMMA2MC total cross section is 39±g pb and the 27NNLO total cross section 
is 44^5 pb, where the uncertainty is dominated by the choice of the nominal scales. 

The comparisons between the experimental cross sections and the predictions by PYTHIA 
and SHERPA are shown in figure 4. In order to compare the shapes of the MC differential 
distributions to the data, their cross sections are rescaled by a factor 1.2 to match the total 
cross section measured in data. PYTHIA misses higher order contributions, as clearly seen 
for low values of A</> 77 , but this is compensated by the parton shower for A</> 77 near n and 
at low Pt,77- It is worth noting that the shoulder expected (and observed) in the PT,-yy 
cross section around the sum of the Ei thresholds of the two photons [34] is almost absent 
in PYTHIA, while SHERPA correctly reproduces the data in this region. This is interpreted 
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as being due to the additional NLO contributions in SHERPA combined with differences in 
the parton showers. Overall, SHERPA reproduces the data rather well, except at large m 77 
and large | cos #* 7 1 . 

The comparisons between the data cross sections and the predictions by 27NNLO 
and DIPHOX+GAMMA2MC are shown in figure 5. In the A<^> 77 ~ ir, low Pt,77 region, 
DIPHOX+GAMMA2MC fails to match the data. This is expected because initial-state soft 
gluon radiation is divergent at NLO, without soft gluon resummation. Everywhere else 
DIPHOX+GAMMA2MC is missing NNLO contributions and clearly underestimates the data. 

With higher order calculations included, 27NNLO is very close to the data within 
the uncertainties. However, the excess at A</> 77 ~ ir and low Pt,^ is still present, as 
expected for a fixed-order calculation. Since the fragmentation component is not calculated 
in 27NNLO, the data is slightly underestimated by 27NNLO in the regions where this 
component is larger: at low A</> 77 , low mass, intermediate Pt,77 (between 20 GeV and 150 
GeV) and large | cos #* 7 1 . 

8 Conclusion 

A measurement of the production cross section of isolated-photon pairs in pp collisions at 
a centre-of-mass energy y/s = 7 TeV is presented. The measurement uses an integrated 
luminosity of 4.9 fb _1 collected by the ATLAS detector at the LHC in 2011. The two 
photons are required to be isolated in the calorimeters, to be in the acceptance of the elec- 
tromagnetic calorimeter (\r/\ < 2.37 with the exclusion of the barrel-endcap transition region 
1.37 < \r}\ < 1.52) and to have an angular separation AR > 0.4 in the 77, (f) plane. Both 
photons have transverse energies E^ > 22 GeV, and at least one of them has > 25 GeV. 

The total cross section within the acceptance is 44.0^42 pb. It is underestimated 
by SHERPA and PYTHIA, which both predict a value of 36 pb with the current ATLAS 
simulation tune. The central value of the cross section predicted by Diphox+gamma2mc, 
39 pb, is lower than the data but it is consistent with data within the theoretical (^g pb) 
and experimental errors. The NNLO calculation of 27NNLO (onnlo = ^-f pb) is m 
excellent agreement with the data. 

The differential cross sections, as a function of the di-photon invariant mass, transverse 
momentum, azimuthal separation and of the cosine of the polar angle of the photon with 
largest transverse energy in the Collins-Soper di-photon rest frame, are also measured. 
Rather good agreement is found with Monte Carlo generators, after rescaling the PYTHIA 
and SHERPA distributions by a factor 1.2 in order to match the integrated cross section 
measured in data and fixed-order calculations, in the regions of phase space studied. All 
generators tend to underestimate the data at large |cos#* 7 |. SHERPA performs rather 
well for most differential spectra, except for high m 77 . PYTHIA is missing higher order 
contributions, but this is compensated by the parton shower for A</> 77 near ir and at low 
Pt,77- I n these same regions the fixed-order calculations do not reproduce the data, due 
to the known infrared divergences from initial-state soft gluon radiation. Everywhere else 
DIPHOX+GAMMA2MC is missing NNLO contributions and clearly underestimates the data. 
On the other hand, with inclusion of NNLO terms, 27NNLO is able to match the data 
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very closely within the uncertainties, except in limited regions where the fragmentation 
component - neglected in the 27NNLO calculation - is still significant after the photon 
isolation requirement. 
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Figure 4. Comparison between the experimental cross sections and the predictions obtained with 
parton-shower LO simulations: m 77 (top left), PT.77 (top right), A</> 77 (bottom left), cos #* 7 (bottom 
right). The LO cross sections have been scaled to the total data cross section, by a factor 1.2. Black 
dots correspond to data with error bars for their total uncertainties, which are dominated by the 
systematic component. The simulated cross sections include only statistical uncertainties. 
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Figure 5. Comparison between the experimental cross sections and the predictions obtained with 
Diphox+GAMMA2mc (NLO) and 27NNLO (NNLO): m 77 (top left), p Tj77 (top right), A<r> 77 (bot- 
tom left), cos (5* 7 (bottom right). Black dots correspond to data with with error bars for their total 
uncertainties, which are dominated by the systematic component. The theoretical uncertainties 
include contributions from the limited size of the simulated sample, from the scale choice and from 
uncertainties on the parton distribution functions and on the hadronization and underlying event 
corrections. 
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A Experimental differential cross section 



The numerical values of the differential cross sections displayed in figures 4 and 5 are quoted 
in tables 2-5. For each bin of the m 77 , Pt,^, A<^> 77 , an d cos#* 7 variables, the cross section 
is given together with its statistical, systematic and total uncertainties. All values are 
divided by the bin width. 





dcr/dm 77 


Statistical 


error 


Systematic 


errors 


Total 


error 


[GeV] 


[pb/GeV] 


high 


low 


high 


low 


high 


low 


[0,20) 


0.0247 


+0.0015 


-0.0015 


+0.0032 


-0.0076 


+0.0036 


-0.0077 


[20,30) 


0.0704 


+0.0032 


-0.0032 


+0.0087 


-0.0140 


+0.0093 


-0.0144 


[30, 40) 


0.091 


+0.004 


-0.004 


+0.011 


-0.015 


+0.012 


-0.015 


[40, 50) 


0.252 


+0.006 


-0.006 


+0.025 


-0.037 


+0.026 


-0.037 


[50,60) 


0.880 


+0.012 


-0.012 


+0.073 


-0.120 


+0.074 


-0.120 


[60, 70) 


0.857 


+0.010 


-0.010 


+0.071 


-0.068 


+0.071 


-0.068 


[70, 80) 


0.626 


+0.008 


-0.008 


+0.051 


-0.052 


+0.051 


-0.053 


[80, 90) 


0.384 


+0.008 


-0.008 


+0.049 


-0.048 


+0.050 


-0.049 


[90, 100) 


0.305 


+0.006 


-0.006 


+0.031 


-0.034 


+0.032 


-0.035 


[100,110) 


0.212 


+0.004 


-0.004 


+0.021 


-0.021 


+0.021 


-0.021 


[110,120) 


0.148 


+0.004 


-0.004 


+0.015 


-0.014 


+0.015 


-0.015 


[120, 130) 


0.122 


+0.003 


-0.003 


+0.015 


-0.010 


+0.015 


-0.011 


[130, 140) 


0.0829 


+0.0025 


-0.0025 


+0.0112 


-0.0073 


+0.0115 


-0.0077 


[140, 150) 


0.0656 


+0.0022 


-0.0022 


+0.0088 


-0.0057 


+0.0091 


-0.0061 


[150, 160) 


0.0535 


+0.0019 


-0.0019 


+0.0072 


-0.0051 


+0.0075 


-0.0054 


[160, 170) 


0.0451 


+0.0017 


-0.0017 


+0.0063 


-0.0038 


+0.0065 


-0.0042 


[170, 180) 


0.0343 


+0.0015 


-0.0015 


+0.0050 


-0.0031 


+0.0052 


-0.0035 


[180, 190) 


0.0262 


+0.0013 


-0.0013 


+0.0032 


-0.0024 


+0.0035 


-0.0027 


[190, 200) 


0.0209 


+0.0011 


-0.0011 


+0.0025 


-0.0019 


+0.0028 


-0.0022 


[200, 225) 


0.0149 


+0.0006 


-0.0006 


+0.0023 


-0.0014 


+0.0024 


-0.0015 


[225,250) 


0.00970 


+0.00049 - 


-0.00049 


+0.00150 - 


-0.00086 


+0.00158 


-0.00099 


[250, 275) 


0.00616 


+0.00039 - 


-0.00039 


+0.00104 - 


-0.00064 


+0.00111 


-0.00075 


[275, 300) 


0.00464 


+0.00036 - 


-0.00036 


+0.00080 - 


-0.00059 


+0.00087 


-0.00069 


[300, 350) 


0.00235 


+0.00017 - 


-0.00017 


+0.00048 - 


-0.00026 


+0.00051 


-0.00031 


[350, 400) 


0.00116 


+0.00011 - 


-0.00011 


+0.00024 - 


-0.00013 


+0.00026 


-0.00017 


[400, 500) 


4.69e-04 


+5.0e-05 - 


-5.0e-05 


+7.9e-05 - 


-4.7e-05 


+9.3e-05 


-6.9e-05 


[500, 800) 


8.6e-05 


+1.3e-05 - 


1.3e-05 


+1.5e-05 - 


-1.0e-05 


+1.9e-05 


-1.6e-05 



Table 2. Experimental cross-section values per bin in pb/GeV for to 77 . The listed total errors 
are the quadratic sum of statistical and systematic uncertainties. 
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£>T,77 


do-/dp Tj77 


Statistical 


error 


Systematic 


errors 


Total 


error 


[GeV] 


[pb/GeV] 


high 


low 


high 


low 


high 


low 


[0,2) 


0.727 


+0.022 


-0.022 


+0.057 


-0.092 


+0.061 


-0.094 


[2,4) 


1.75 


+0.04 


-0.04 


+0.13 


-0.23 


+0.13 


-0.23 


[4,6) 


2.03 


+0.04 


-0.04 


+0.15 


-0.23 


+0.15 


-0.23 


[6,8) 


1.88 


+0.04 


-0.04 


+0.15 


-0.21 


+0.16 


-0.21 


[8, 10) 


1.72 


+0.03 


-0.03 


+0.13 


-0.19 


+0.14 


-0.19 


[10,12) 


1.40 


+0.03 


-0.03 


+0.12 


-0.16 


+0.12 


-0.16 


[12, 14) 


1.28 


+0.03 


-0.03 


+0.10 


-0.13 


+0.11 


-0.13 


[14, 16) 


1.122 


+0.026 


-0.026 


+0.093 


-0.114 


+0.097 


-0.117 


[16,18) 


0.999 


+0.024 


-0.024 


+0.086 


-0.090 


+0.090 


-0.093 


[18,20) 


0.810 


+0.021 


-0.021 


+0.072 


-0.076 


+0.075 


-0.079 


[20,25) 


0.674 


+0.012 


-0.012 


+0.056 


-0.074 


+0.058 


-0.075 


[25,30) 


0.492 


+0.011 


-0.011 


+0.041 


-0.045 


+0.043 


-0.047 


[30,35) 


0.405 


+0.009 


-0.009 


+0.034 


-0.043 


+0.035 


-0.044 


[35, 40) 


0.325 


+0.009 


-0.009 


+0.028 


-0.034 


+0.030 


-0.035 


[40, 45) 


0.272 


+0.008 


-0.008 


+0.024 


-0.027 


+0.026 


-0.028 


[45, 50) 


0.282 


+0.008 


-0.008 


+0.023 


-0.027 


+0.024 


-0.028 


[50,55) 


0.235 


+0.007 


-0.007 


+0.023 


-0.025 


+0.025 


-0.026 


[55,60) 


0.194 


+0.006 


-0.006 


+0.019 


-0.024 


+0.021 


-0.024 


[60,65) 


0.150 


+0.006 


-0.006 


+0.015 


-0.016 


+0.016 


-0.017 


[65, 70) 


0.102 


+0.005 


-0.005 


+0.013 


-0.012 


+0.014 


-0.013 


[70, 75) 


0.0836 


+0.0041 


-0.0041 


+0.0103 


-0.0087 


+0.0111 


-0.0096 


[75,80) 


0.0748 


+0.0036 


-0.0036 


+0.0087 


-0.0086 


+0.0094 


-0.0093 


[80, 90) 


0.0521 


+0.0021 


-0.0021 


+0.0059 


-0.0056 


+0.0063 


-0.0059 


[90, 100) 


0.0381 


+0.0017 


-0.0017 


+0.0043 


-0.0036 


+0.0047 


-0.0040 


[100,110) 


0.0239 


+0.0013 


-0.0013 


+0.0028 


-0.0023 


+0.0031 


-0.0026 


[110,120) 


0.0175 


+0.0011 


-0.0011 


+0.0024 


-0.0016 


+0.0027 


-0.0019 


[120, 130) 


0.0106 


+0.0009 


-0.0009 


+0.0015 


-0.0011 


+0.0017 


-0.0014 


[130, 140) 


0.0090 


+0.0008 


-0.0008 


+0.0012 


-0.0008 


+0.0015 


-0.0012 


[140, 150) 


0.00646 


+0.00064 - 


-0.00064 


+0.00089 - 


-0.00063 


+0.00110 


-0.00090 


[150, 175) 


0.00333 


+0.00031 - 


-0.00031 


+0.00047 - 


-0.00039 


+0.00056 


-0.00049 


[175, 200) 


0.00195 


+0.00023 - 


-0.00023 


+0.00025 - 


-0.00017 


+0.00034 


-0.00028 


[200, 250) 


0.00077 


+0.00010 - 


-0.00010 


+0.00012 - 


-0.00008 


+0.00016 


-0.00013 


[250, 500) 


1.18e-04 


+1.7e-05 - 


1.7e-05 


+1.8e-05 - 


-1.2e-05 


+2.5e-05 


-2.1e-05 



Table 3. Experimental cross-section values per bin in pb/GeV for pt^ 11 . The listed total errors 
are the quadratic sum of statistical and systematic uncertainties. 
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A0 77 


dc/dA^p, 


Statistical 


error 


Systematic 


errors 


Total 


error 


[rad] 




[pb /rad] 


high 


low 


hieh 

o 


low 


high 


low 


[0.00, 


0.50) 


2.68 


+0.08 


-0.08 


+0.22 


-0.52 


+0.24 


-0.52 


[0.50, 


1.00) 


3.10 


+0.09 


-0.09 


+0.25 


-0.36 


+0.26 


-0.37 


[1.00, 


1.50) 


3.46 


+0.09 


-0.09 


+0.29 


-0.36 


+0.31 


-0.37 


[1.50, 


1.75) 


4.51 


+0.15 


-0.15 


+0.36 


-0.42 


+0.39 


-0.44 


[1.75, 


2.00) 


6.26 


+0.17 


-0.17 


+0.53 


-0.59 


+0.55 


-0.61 


[2.00, 


2.25) 


8.93 


+0.20 


-0.20 


+0.73 


-1.02 


+0.76 


-1.04 


[2.25, 


2.35) 


11.6 


+0.4 


-0.4 


+0.9 


-1.1 


+1.0 


-1.2 


[2.35, 


2.45) 


13.9 


+0.4 


-0.4 


+1.1 


-1.3 


+1.1 


-1.4 


[2.45, 


2.55) 


17.4 


+0.4 


-0.4 


+1.3 


-1.6 


+1.4 


-1.7 


[2.55, 


2.65) 


21.8 


+0.5 


-0.5 


+2.0 


-2.2 


+2.0 


-2.3 


[2.65, 


2.70) 


26.7 


+0.8 


-0.8 


+2.1 


-2.4 


+2.3 


-2.5 


[2.70, 


2.75) 


29.6 


+0.8 


-0.8 


+3.2 


-3.4 


+3.3 


-3.5 


[2.75, 


2.80) 


35.8 


+0.9 


-0.9 


+2.8 


-3.1 


+3.0 


-3.3 


[2.80, 


2.85) 


42.9 


+1.0 


-1.0 


+3.3 


-3.8 


+3.4 


-3.9 


[2.85, 


2.90) 


48.4 


+1.1 


-1.1 


+3.8 


-4.4 


+3.9 


-4.6 


[2.90, 


2.95) 


57.4 


+1.2 


-1.2 


+4.3 


-4.6 


+4.4 


-4.7 


[2.95, 


3.00) 


71.7 


+1.3 


-1.3 


+5.4 


-5.9 


+5.6 


-6.1 


[3.00, 


3.05) 


80.8 


+1.4 


-1.4 


+6.1 


-7.4 


+6.3 


-7.6 


[3.05, 


3.10) 


100.5 


+1.6 


-1.6 


+7.3 


-8.5 


+7.4 


-8.6 


[3.10, 


3.14) 


107.6 


+1.8 


-1.8 


+7.9 


-9.6 


+8.1 


-9.8 



Table 4. Experimental cross-section values per bin in pb/rad for A0 77 . The listed total errors 
are the quadratic sum of statistical and systematic uncertainties. 
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cos & 77 


dcj/acos x) 


statistical 


error 


Q„ .-.4- „™ „ j.: ., 

systematic 


errors 


lotal 


error 




[P b J 


high 


low 


high 


low 


high 


low 


[—1.00, —0.92) 


11.0 


+0.4 


—0.4 


+ 1.2 


— 1.6 


+ 1.3 


— 1.7 


[—0.92, —0.84) 


1 A f 

14.5 


+0.4 


—0.4 


+ 1.1 


— 1.6 


+ 1.2 


— 1.6 


[—0.84, —0.76) 


16.5 


+0.5 


—0.5 


+ 1.3 


— 1.6 


i 1 A 

+ 1.4 


— 1.6 


[—0.76, —0.68) 


18.8 


+0.5 


—0.5 


i 1 A 

+1.4 


— 1.6 


+1.5 


— 1.7 


[—0.68, —0.60) 


21.8 


+0.6 


—0.6 


+1.5 


— 1.9 


+1.6 


— 1.9 


[—0.60, —0.52) 


22.7 


+0.6 


—0.6 


+1.8 


—2.0 


+1.9 


—2.1 


[—0.52, —0.44) 


23.9 


+0.6 


—0.6 


+1.8 


O 1 

—2.1 


+ 1.9 


O 1 

—2.1 


[—0.44, —0.36) 


24.6 


+0.6 


—0.6 


+1.8 


—2.0 


+1.9 


— 2.1 


[—0.36, —0.28) 


25.2 


+0.6 


—0.6 


+1.8 


—2.2 


+1.9 


—2.3 


[—0.28, —0.20) 


28.3 


+0.6 


—0.6 


i O 1 

+2.1 


O f* 

—2.6 


i O O 

+2.2 


O T 

—2.7 


[—0.20, -0.12) 


28.1 


+0.7 


—0.7 


+2.0 


—2.6 


+2.1 


—2.6 


[—0.12, —0.04) 


29.6 


+0.7 


—0.7 


+2.3 


O T 

—2.7 


+2.4 


—2.8 


r r\ r\ A c\ r\ a\ 

[—0.04, 0.04) 


31.4 


+0.7 


—0.7 


+2.5 


—2.9 


+2.6 


—3.0 


[0.04, 0.12) 


29.0 


+0.7 


—0.7 


+2.4 


o r 

—2.5 


+2.5 


—2.6 


[U.lz, U.zU J 


Z 1 .0 


+U. i 


n v 
— U. / 


l 1 
+ Z.1 


— o.U 


1 o o 


— 6.1 


[0.20,0.28) 


26.2 


+0.6 


-0.6 


+1.8 


-2.3 


+2.0 


-2.3 


[0.28,0.36) 


25.7 


+0.6 


-0.6 


+2.0 


-2.1 


+2.1 


-2.2 


[0.36,0.44) 


25.0 


+0.6 


-0.6 


+1.9 


-1.9 


+2.0 


-2.0 


[0.44,0.52) 


24.9 


+0.6 


-0.6 


+1.7 


-2.0 


+1.8 


-2.1 


[0.52,0.60) 


22.7 


+0.6 


-0.6 


+1.9 


-2.2 


+2.0 


-2.2 


[0.60,0.68) 


21.3 


+0.6 


-0.6 


+1.6 


-1.8 


+1.7 


-1.9 


[0.68,0.76) 


19.5 


+0.5 


-0.5 


+1.4 


-1.8 


+1.5 


-1.9 


[0.76, 0.84) 


16.9 


+0.5 


-0.5 


+1.4 


-1.6 


+1.5 


-1.7 


[0.84, 0.92) 


13.9 


+0.4 


-0.4 


+1.1 


-1.7 


+1.2 


-1.8 


[0.92, 1.00) 


11.7 


+0.4 


-0.4 


+1.0 


-1.4 


+1.1 


-1.4 



Table 5. Experimental cross-section values per bin in pb for cos#* 7 . The listed total errors are 
the quadratic sum of statistical and systematic uncertainties. 
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